Simulating Wess-Zumino Supersymmetry Model in Optical Lattices 
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We study a cold atom-molecule mixture in two-dimensional optical lattices. We show that by 
fine-tuning the atomic and molecular interactions, Wess-Zumino supersymmetry (SUSY) model in 
2-|-l-dimensions emerges in the low-energy limit and can be simulated in such mixtures. At zero 
temperature, SUSY is not spontaneously broken, which implies identical relativistic dispersions of 
the atom and its superpartner, bosonic diatom molecule. This defining signature of SUSY can 
be probed by single-particle spectroscopies. Thermal breaking of SUSY at a finite temperature is 
accompanied by a thermal Goldstone fermion, i.e., phonino excitation. This and other signatures 
of broken SUSY can also be probed experimentally. 



PACS numbers: 67.85. Pq, 37.10.Jk, ll.30.Pb 



Introduction. - Wess and Zumino proposed the first 
space-time supersymmetry (SUSY) model (WZ-SUSY 
model) 36 years ago [1|. Since then SUSY has become 
a fundamental ingredient of theories beyond the stan- 
dard model in high-energy physics [2|. However, none of 
super partners of the known elementary particles have 
been found thus far; it remains to be seen if they can be 
detected in the energy range of Large Hadron Collider. 

On a different front, nonrelativistic SUSY (a Bose- 
Fermi symmetry unrelated to space-time symmetry) has 
attracted considerable recent interest in the cold atom 
community, as it can be realized by using Bose-Fermi 
atom (molecule) mixtures which are loaded in optical lat- 
tices. Examples include attempts to simulate the nonrel- 
ativistic limit of superstring by trapping fermionic atoms 
in the core of vortices in a Bose-Einstein condensate [3| ; 
study of the SUSY effect in an exactly solvable one- 
dimensional Bose-Fermi mixture with Bethe ansatz [4J; 
and SUSY models for nonrelativistic particles in various 
dimensions 5->7] . In Ref . [6] , we studied perhaps the sim- 
plest cold atom SUSY model and discussed detecting the 
Goldstino-like mode due to SUSY breaking by measur- 
ing a single fermion spectral function. In a further work 
[7|, we developed a SUSY response theory to photoasso- 
ciation in a cold fermionic atom system. Although these 
studies are interesting and some results may be general 
for many SUSY systems [7] , SUSY in these nonrelativistic 
systems is very different from the relativistic (or space- 
time) SUSY in high-energy physics. 

In this Letter, we propose a way to simulate the sim- 
plest relativistic SUSY model, the WZ-SUSY model Q. 
We show that it can emerge in the low-energy limit of 
a cold atom-molecule mixture in properly chosen two- 
dimensional lattices. The first requirement is the exis- 
tence of Dirac points in the Brillouin zone. Recently, 
such models based on a honeycomb lattice or graphene- 
like structure have been proposed [8|. In these models, 
two Dirac points K and K' are related to each other 
by K' = —K, as required by time-reversal symmetry. 
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This means This means that two fermionic atoms which 
form a usual BCS pair or a diatom molecule belong to 
two different Dirac points. To simulate the WZSUSY 
model, however, one needs a Klein-Gordon field as the 
Dirac fermions superpartner which corresponds to a di- 
atom molecule made by two Dirac fermions from the 
same Dirac point. Such molecules carry a 2K ^ mo- 
mentum and are energetically unfavorable as a result. In 
a recent work, Lee attempted to avoid this difficulty by 
introducing frustrated hopping for the molecules, such 
that the boson dispersion has minima at ±2X instead of 
zero [9|. It is found that the massless WZ-SUSY model 
emerges at the boson's superfluid-insulator critical point. 

In this work, we show that the WZ-SUSY model can 
emerge not only at the critical point. We use a lattice 
model studied recently by Liu et al. 10] instead. This 
is a square lattice model in which the Dirac points at 
K — (0,0) and K' = (0,7r) are their own negatives, 
as (0, — tt) = (0,7r). This means a diatom molecule 
made of two atoms from the same Dirac points has zero- 
momentum for 2K = (0,0) and 2K' = (0, 27r) ee (0,0). 
With this setup, we can simulate the WZ-SUSY model 
more straightforwardly, after appropriate interactions are 
introduced and fine-tuned. 

The research interest in WZ-SUSY models has been 
renewed recently [ll|. No spontaneous breaking of the 
SUSY implies there are equal poles in the single-particle 
spectral functions of both the Dirac field and the Klein- 
Gordon field. A further calculation showed that these 
single-particle spectral functions are not renormalized 
from their free particle ones 
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of the SUSY and may be detected by the established 
techniques of the single-particle spectroscopies [13J- It 



is known that a thermal bath always breaks the SUSY 
14| and this thermal breaking of the SUSY is accompa- 
nied by a thermal Goldstone fermion, phonino 12]; thus 
studying this model at finite temperature sheds light on 
physics of SUSY breaking. 

There are many studies of SUSY in space-time lat- 



tice models [15|. The significant difference between the 
the present work and those lattice SUSY models is that 
while the latter are supersymmetric on the lattices, we 
study the emergence of SUSY from a microscopic space 
lattice (but continuous time) model with no SUSY to 
begin with. 

Free Fermion Lattice Model and Continuum Limit. - 
We briefly recall the lattice model proposed in Ref. |lO|. 
Consider a single-component fermionic atom gas loaded 
in a square lattice. The potential minimum in the sublat- 
tice A is higher than that in the sublattice B. Two states 
with the energy difference 2M , the s-orbital at the ^-sites 
and the p-orbital at the i?-sites, form a pseudospin-1/2 
subspace. The sublattices are anisotropic with 1,2,3, and 
4, the next nearest neighbor sites (Fig. 1(b)). The hop- 
pings between the nearest and next nearest sites are taken 
into account. The corresponding hopping amplitudes are 

tA,A+5^(y^ = ^tA,A-S^(y^ = tAB, ^Al = tA3, tA2 — tA4, 

tsi = ifis, and tB2 = tsA with 5x ( 5y) being the unit 
vector in the x- (y-) direction. In addition, a periodic 
gauge field generated by two opposite-traveling standing 
wave laser beams coupling with atoms 16] is introduced. 
This gives rise to a tunable staggered Peierls phase ±0q 
along the vertical links and vanishing in the horizontal 
and 1, 2, 3, 4 links. With these the single-fermion Hamil- 
tonian is given by 



i/(k) = PxO^)(yx +Py{\i.)(Ty + h^{\<i)a^ 



(1) 



where p^ = 2t^s sin0osin(/cj^a), py =^ 2i^B(sin(A:2;a) + 
cos^o sin(fcya)) and h^ — —M — tg cos(A;j.a) cos(fcj,a) — 
2f sin(fc3,a) sin(fcj^a) with to = tj^i — t^i + tA2 ^ ^B2 and 
i = {tAi - tsi + ts2 - iA2)/2. When M = ±to 7^ 0, 
there is a unique gapless Dirac point: either K = (0, 0) 
or K' = (0,7r). We choose Oo — 7r/2 and define the 
'speed of light' Vs = 2tABCi. In the continuum limit and 
near the Dirac points, Px{K + 5k) ^ 2tABCiSky = VsQx 
and Px{K' + Sk) ^ —2tAB0-Sky — —VsQx ', Py{K + 6k) ^ 
2tABaSkx = VsQy and Py{K' + Sk) ~ 2tABaSkx = v^qy', 
and h^{K + 6k) = -M - to ^ tuq and h^{K' + 6k) = 
—M + to = "^TT- Thus, for the Dirac fields ^(x) near K 
and C(a;) near K' , the effective Hamiltonian reads 



i/(0) = V, d^xCH-^<da + moa,)^ 



+ Vs d^x(^{-ia"ida + mT,az)C„ 



(2) 



t^ — ±0-^ and a\ 
field ^(r) is given by 



a^; the two-component Dirac 
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with 6c, 



^q - ^^^1^ = -0-q and similarly for C- The 
momentum cut-off A corresponds to that when lattice 
fermion dispersion deviates severely from the linear one. 
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FIG. 1: (a) Josephson tunneling between the atom-molecule 
mixture (lower lattice plane) and the dimolecule Bose- 
Einstein condensate nearby (upper plane). The orange dots 
are molecules in the mixture and red dots are dimolecules. 
Fermionic atoms are in the lattice sites, (b) The square lat- 
tice where 1,2,3,4 denote the next nearest neighbor sites. 



The mass terms here can be fine-tuned. When M = ±io 
which is not zero in this lattice setup [Ifl], one of the 
Dirac field is massless and another is massive. The lat- 
ter can be integrated out in the low-energy limit. We 
note the zero matter density (necessary for Lorentz in- 
variance) in relativistic quantum field theory corresponds 
to fermionic atoms being at half filling in this lattice re- 
alization. After an external source is introduced, fermion 
number (including those forming molecules) will fluc- 
tuate but average at half filling. To facilitate pairing 
or molecule formation, we introduce attraction between 
fermionic atoms, which is modeled in a two-channel fash- 
ion below. For such spinless fermions, the two-atom at- 
traction and p-wave type bound state have already been 
achieved experimentally 17| . 

Two-channel Model. - We take toq = and ttt.^ ^ 0, 
and integrate out C. in the low-energy limit where only the 
states with their energy lower than min{TO7r, E\} are rel- 
evant. We now extend our Hamiltonian to a two-channel 
model, i.e., the lowest two hyperfine atom states with 
two-atom scattering states in open channel and the two- 
atom bound state (Feshbach molecule) in closed channel. 
We denote C*^°-'(a;) the Dirac fermions in the open chan- 
nel and C (2^) the Dirac fermion in the closed channel. 
Analogous to the many body theory of the atom- molecule 
coherence in Ref. Ig, the effective Lagrangian describing 
this two-channel Dirac fermion model is given by 
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'T^s^.e 
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_Jj(c)£(c)t£(c)t£(c)£(c) 






where ct^ = {L,crx,<Ty) and d^ — (9t,UsV). U^'^^ and 
jjico) g^j.g ^Yic interaction between closed channel fermions 
and the interchannel interaction, respectively. We have 
neglected the background interaction in open channel. 
By introducing the pairing field A(r, t) for ^(''^ via 
a Hubbard-Stratonovich transformation and integrating 
out ^('^', the resulting Lagrangian is given by 



/:(^(°),A) = -ic(°)t^j.^Ma^^(o) 






TrlnG(")-i. 



(3) 



The inverse of the propagator G^'^^ of ^^'^^ is given by 



q(c)-1 ^ 



^o^'^^ 



J^d., 
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St 



where S 



.^(co)^(o)^(o)^ 



Expanding the Lagrangian 



in powers of A and its gradients yields 
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df.^^d^'ip 
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l"^!' 



^i'+o(i^n, 



(4) 

where if oc A/U^'^^ is the Feshbach molecular field with 
the detuning energy £„ and the interacting strength 
A ex (C/''^))^. We have Vb = Vs in the weak coupHng 
limit (i.e., f/^'^'™^ much smaller all other energy scales 
in the system including rrij^ and E\) due to (emergent) 
Lorentz invariance. Lattice effects (which break Lorentz 
invariance) give rise to nonuniversal corrections to Vg and 
Vb] thus tuning of one parameter {e.g., molecule disper- 
sion through an additional lattice potential seen by the 
molecule only) is needed to ensure Vb = Vs to maintain 
Lorentz invariance in the low-energy limit. Also included 
in ^ is the Yukawa coupling between (p and ^(°\ i.e., 

^^C = -f (^e^°^^^l°^^ + ^^d°^e^°^) with g (X -2C/(-). 

WZ-SUSY Model: Massless. - For simplicity, we drop 
the superscript of ^^°'^ hereafter. By combining ([3]), 
((4]) and the Yukawa coupling together, the effective La- 
grangian after neglecting 0(|</3|^) is given by 



>c(e,^) 



-2^m'/'^^^<^- 2^ 



^1 



^|2 - *eV''5^e 

'-kv^kl + v^^ii2). (5) 



Tuning Em = by varying U^'^' , and further tuning 
pair-pair (or molecule-molecule) interaction by varying 
[7*^'^°-' so that the coupling constant \ — g^, the effec- 
tive Lagrangian £(^, tp) is exactly the massless WZ-SUSY 
model with the SUSY under the SUSY transformations 
5ip — eV^<^, and 55, = a^^a^ed^^p'^ — f (yS^^e where e is a 
constant two-component spinor parameter. 

WZ-SUSY Model: Massive. - To have a massive WZ- 
SUSY model, we need to introduce an external source. 
This can be realized by putting a Bose-Einstein con- 
densate of dimolecules nearby, which is made of pairs 
of molecules (or 4-atom molecules) (see Fig. 1(a)). 
Through Josephson tunneling with an amplitude n, the 
dimolecule condensate exchanges pairs of molecules with 
the mixture. The effective Lagrangian reads 

c{i,ip,^) = -\d,p:^d^p>^i5^a^d^i 






(6) 



where '^ is the external dimolecular field. There is a 
global U(l) symmetry (called R-symmetry) under ^ — >■ 



e^^i, ip -> e^'^^ip and * -> e'*''^* [ll|. If * slowly varies 



in space-time, it is also SUSY invariant under 5if = e^cr^^ 



and 5S_ = a^a'^ed^i.tf'^ 



f^t2. 



4k*' 



e. By taking $ to 



be its condensed order parameter (^) = (^pt) = m^/8K, 
the R-symmetry is broken and reduced to a discrete Z2 
symmetry with f — >■ i^ and (/?—>■ —t^, and the on-shell 
WZ Lagrangian appears (up to an additive constant) 



C{^,(p,Tn) 



Id^^^d'^p-i^^a'^d^^ (7) 



5 / +2 '^^\f -2 



^2 



The SUSY is exact by replacing vl/t with (*t) in the 
SUSY transformations. The Z2 symmetry is always spon- 
taneously broken in one of the degenerate ground states 
with ip = (j) ± in/ g. The SUSY Lagrangian with sponta- 
neous breaking of Z2 becomes 



L = - 



*^V^a^eT^™(^kI + ei6 



9 



('/'Ckl + 0ki6)- 



^101^ 



(8) 



This is the 2+1-dimensional reduction of the original 
WZ-SUSY model in 3-^1 dimensions [l|,l^. 

Supercurrent and Supercharge. - SUSY leads to a 
conserved supercurrent, whose conserved supercharges 
are generators of SUSY. The supercurrent is defined by 
5 j dtd^xC = J dtd^xe^a'^df^Jj^. The supercharges are 
then given by Q = Jd'^xjf{x) and Qt. The SUSY 
transformation generated by Q for a field O reads SO = 
—ie^(jy[Q, 0]±. We focus on the on-shell model with the 
Z2 symmetry spontaneously broken, where the on-shell 
supercurrent [19| is given by J^ = ia'^a'^^d^4> + *§ (^^ -'- 
2m(t)/ g)ay(j^^. The SUSY spontaneous breaking is sig- 
naled by {{Q, O}) ^ for a fcrmionic operator O. How- 
ever, for this simplest SUSY model, the SUSY is not 
spontaneously broken at zero temperature [2|. 

Nonrenormalization . - In this simplest WZ-SUSY 
model, single-particle Green's functions are not renor- 
malized due to (unbroken) SUSY. For example, the 
renormalization to Klein-Gordon field's propagator in a 
one-loop self-energy calculation is given by q^ — m^ — > 
g2 - mliq) with [12] m^q) ^m + gR{A)o + Oigl{A)l) 
where A — Ke (j>. For 2+1 dimensions, due to the 
nonzero anomalous critical exponents [20[, the coupling 
constant may be renormalized to gu. However, {A)q cx 
{{QiS.}}o — because the SUSY is not spontaneously 
broken. The mass of the Dirac field is also not renor- 
malized as required by SUSY. Therefore, the single par- 
ticle Green's functions, both of the Dirac and Klein- 
Gordon fields, are not renormalized from their free ver- 
sion. The spectral functions of the Green's functions can 



be measured by the single-particle spectroscopic tech- 
nique which has been developed recently (131] . The non- 
renormalization of the Green's functions implies sharp 
peaks in their spectral functions, with identical relativis- 
tic dispersions for the atoms and molecules. Experimen- 
tally this would be the hallmark of achieving SUSY. 

Thermal breaking of SUSY. - By replacing t by zt, 
the imaginary time, the Euclidean version of Lagrangian 
(HI) describes WZ-SUSY model in finite temperature 
T. When T ^ 0, SUSY is always broken because 
{{Q, Q^ay})T = {'^''PflT ^ with P^ being the energy- 
momentum operator |14| . due to the nonvanishing ther- 
mal energy. This SUSY thermal breaking is accompanied 
by a thermal Goldstone fermion (phonino) but not nec- 
essarily by a phonon because the Lorentz symmetry is 
also broken by (Po)t t^ 0[12]. The phonino dispersion is 
given by[l4| qo = ±Vss\(l\ where the SUSY sound velocity 
Vss = Vs/B for T ^ m and Vss — Tvs/m for T <^m. 

To detect the phonino mode, one can consider the 
response to an external 'fermionic' field coupled to 
the supercurrent. The phonino is a pole of the 
supercurrent-supercurrent correlation function. This ex- 
ternal 'fermionic' field can be a combination of an exter- 
nal photon with another hyperfine state of the fermionic 
atom which is decoupled to the mixture. We have studied 
this kind of SUSY response theory for a nonrelativistic 
SUSY mixture [TJ. However, the difficulty in the present 
case is that the supercurrent is not so simple as that in 
the nonrelativistic theory, and thus the coupling between 
the external 'fermionic' field and the supercurrent is not 
that easy to be experimentally handled. 

Replacing {A)q by (^)t, the masses are thermally 
renormalized. The masses of A, i? (0 = yl + iB) and 
the spinor ^ have been calculated in low temperature and 
high temperature limits [12|. Namely, in 2+1-dimensions 
up to one- loop, for T <^m, one has ms = ?7i, m\ — m?g ex 

where a = 2!I-g^-™-/T. f^j. 

2g2T, ml=m'^ - g^T. 
These unequal masses of these fields signal SUSY break- 
ing, and can be probed quantitatively. In particular, 
we spectroscopy measurements to show double peaks in 
the molecule spectral function due to the unequal masses 
between A and B components, while the atom spectral 
function has a single peak with a mass of the Dirac field 
equal to neither ttla nor niB- 

Experimental Challenges. - Optical lattices that trap 
cold atoms can be routinely set up in laboratory. The 
staggered Peierls phase originates from production of 
the artificial magnetic field [l6(. As discussed earlier, 
one needs to tune three parameters to achieve SUSY 



g'^ma, m^ — wla oc g^ma 
T ^ m, TUB = rn, m\ — rn^ 



21| : atom-atom interaction, molecule-molecule interac- 
tion, and molecule velocity. The former two may be 
done by adjusting the real magnetic field in Feshbach 
resonance while the latter may be related to the inter- 
action between laser field and the Feshbach molecule. 
These are all achievable within existing experimental ca- 



pabilities. Perhaps the biggest challenge is finding the 
right fermionic atom, which needs to have a highly tun- 
able interaction through a p-wave Feshbach resonance. 
It also needs to support a sufficiently stable dimolecule 
(or 4-atom bound) state, whose condensate provides the 
source term in Eq. ([5]), which gives rise to equal parti- 
cle masses. Experimentally, one needs to overcome the 
atom loss due to the heating of the atom gas caused by 
three- and four-body collisions. Without the last ingre- 
dient however, one can still realize the massless version 
of the WZ-SUSY model, Eq. ©(with e^ tuned to zero), 
which already contains very rich SUSY physics. De- 
spite these and other challenges, we believe simulating 
the WZ-SUSY model using cold atom-molecule mixtures 
is a worthwhile endeavor, as as generalizing the present 
model to a (3-|-l)-dimensional model is straightforward, 
and then it provides new opportunities to explore the real 
space-time SUSY physics. 
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